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Abstract 

Aims. We observed the r Boo system with the HARPS-N spectrograph to test a new observational strategy aimed at jointly studying 
asteroseismology, the planetary orbit, and star-planet magnetic interaction. 

Methods. We collected high-cadence observations on 11 nearly consecutive nights and for each night averaged the raw FITS files 
using a dedicated software. In this way we obtained spectra with a high signal-to-noise ratio, used to study the variation of the Ca n 
H&K lines and to have radial velocity values free from stellar oscillations, without losing the oscillations information. We developed 
a dedicated software to build a new custom mask that we used to refine the radial velocity determination with the HARPS-N pipeline 
and perform the spectroscopic analysis. 

Results. We updated the planetary ephemeris and showed the acceleration caused by the stellar binary companion. Our results on 
the stellar activity variation suggest the presence of a high-latitude plage during the time span of our observations. The correlation 
between the chromospheric activity and the planetary orbital phase remains unclear. Solar-like oscillations are detected in the radial 
velocity time series: we estimated asteroseismic quantities and found that they agree well with theoretical predictions. Our stellar 
model yields an age of 0.9 ± 0.5 Gyr for r Boo and further constrains the value of the stellar mass to 1.38 ± 0.05 M 0 . 
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* Based on observations made with the Italian Telescopio Nazionale 
Galileo (TNG) operated on the island of La Palma by the Fundacion 
Galileo Galilei of the INAF at the Spanish Observatorio Roque de los 
Muchachos of the IAC in the frame of the program Global Architecture 
of the Planetary Systems (GAPS). 

** Table 1 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
http://cdsweb.u-strasbg.fr/cgi-bin/qcat7J/A-l-A/ 


1. Introduction 

While the number of confirmed exoplanets is rapidly increasing, 
we are improving our ability to characterize them by studying 
their parent stars and the star-planet tidal or magnetic interac¬ 
tions. At the moment, the most interesting targets for characteri¬ 
zation are transiting exoplanets, but the large population of non¬ 
transiting planets solicits new methods and observation strate¬ 
gies for this purpose. 

With the intent of characterizing planetary systems with 
a spectro-photometric approach, we selected the well-known 
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system r Bootis A (HD 120136, F 6 V, V=4.49) as a test 
case in the Global Architecture of Planetary System^] (GAPS, 
ICovino et al .1120131 ) programme. GAPS is a joint effort of Italian 
researchers in collaboration with a few experts abroad, and the 
GAPS team manages a long-term observation al program with 
the h igh-precision HARPS-N spectrograplfl ( Cosentino et al. I 
120121) at the Telescopio Nazionale Galileo (TNG), r Boo A is 
included in the subprogram dedicated to the characterization of 
planetary systems through studies of the interactions between 
the planets and their central stars, r Boo’s brightness allows for 
high-resolution spectroscopy and asteroseismology with a lim¬ 
ited investment of telescope time. 

The bright F 6 V star has a faint M2V companion (r Boo B, 
separation 1.83”. fDrummondll2014l) . forming a long-period bi¬ 
nary system. The primary component (hereafter r Boo for the 
sake of uniformity with current literature) was firs_t_claimed to 
host a planet with a period of 3.312 days bv lButler et al.l dl997h : 
they used the radial velocity (RV) method. Thanks to its bright¬ 
ness, this system was employed in the past to develop new analy¬ 
sis techniques that have later been applied to other stars. Among 
the most important discoveries are the definition of upper lim¬ 
its on reflected starlight,_which provid es a maximum v alue for 
the planet’s albedo (iCollier Cameron et al.l 119991 iRodler et al.l 
120101) . the spectroscopic detection of CO absorption lines in 
the planet atmosphere, which permitted determining the incli¬ 
nation angle of the system and thus the exact mass of the plane t 
(M p =5.95 ± 0.28 M Jup , iBrogi et al.ll2()H iRodler et al.lDOTl) 
and, very recently, the first detection of water vapor in the at mo- 
sphere of a non-transiting exoplanet dLockwood et alJl2014lf . 

r Boo was also considered in searches for effects of star- 
planet magnetic interaction (SPMI). A peculiar characteristic of 
this star is the optically variable modulation of the light curve 
that is probably due to photospheric spots that persisted at fixed 
longitudes for_a_few hundred days, as observed by the MOST 
satellite (IWalker et al.120081) . SPMI is a long-debated issue, with 
the best evidence coming from a modulation of chromospheric 
activity tracers phased with the pla netary orbital period rather 
than with the stellar rotation period (IShkolnik et al.l[2005l l2008fc 
lTanzall2009ll20T2lf . 

Assessment of this behavior requires long-term monitoring 
of stars with close-in massive planets (hot Jupiters), very high 
signal-to-noise spectra (S/N> 300 at 3950A), and adequate re¬ 
solving power (R > 80000) to measure variability in the core 
of deep chromospheric lines such as the Can H&K doublet. 
Detecting periodicities equal to the planetary period is crucial 
because variability can be due to intrinsic stellar activity, not re¬ 
lated to the presence of hot Jupiters. On the other hand, the study 
of stellar activity is important as a source of noise in the search 
and characterization of new planets in extra-solar systems. 

In the case of r Boo, pre vious sear ches for SPMI in the Ca n 
H&K lines were ambiguous (IShkolnik et al .12005! ). One possible 
reason is that SPMI is due to magnetic stresses between the stel¬ 
lar and planetary magnetic fields, but this effect is very limited in 
r Boo because the stellar rotation is known to b e sync hronized 
with the orbital motion of the planet (e.g., IHenrv et al.l!2000l) . 
Nonetheless, this system remains an interesting target because 
the parent star shows ev idence of magnetic cycles, with yearly 
polarity reversal s (e.g., iGatala et aD 20071 : iDonati et~ai]l2008t 
iFares et akll2009l) . 


1 http://www.oact.inaf.it/exoit/EXO-IT/Projects/Entries/2011/12/ 
27_GAPS.html 

2 http://www.tng.iac.es/instruments/harps/ 


Asteroseismology of stars hosting exoplanets received a 
great boost from the photometric time series collected by space 
missions. On the other hand, spectroscopic campaigns aimed at 
detecting solar-like oscillations are very hard to organize be¬ 
cause a long time baseline is required Jo resolve and identify 
the excited modes. For instance. iBazot et al.l (1201 2l) were able to 
obtain a measure of the large separation of 18 Sco by means 
of 2833 data points collected over 12 complete nights. It is 
difficult to insert an asteroseismic program in a large project 
such as GAPS, which covers a multiplicity of subprograms and 
goals with limited telescope time. Moreover, the exposure times 
have to be very short to monitor the solar-like oscillations, and 
consequently, the asteroseismic targets have to be bright stars. 
Therefore we used the scientific case of r Boo as a pathfinder, ap¬ 
plying the high-cadence observational strategy for asteroseismic 
purposes and matching it with some exoplanetary goals (search 
for additional planets and study of the star-planet interaction). 
The paper is structured as follows: Sect. [2] presents the obser¬ 
vations and describes the data reduction. Section [3] is devoted 
to the spectral analysis, with results presented in Sect. [4] (stellar 
parameters). Sect. 0 (orbital parameters). Sect. [ 6 ] (stellar activ¬ 
ity), Sect.Rasteroseismology), and Sect. Revolutionary stage). 
Conclusions are presented in Sect. [9] 


2. Observations and data reduction 

t Boo was observed with HARPS-N on 11 nights between April 
13 and May 8 , 2013. Very good phase coverage of the orbital 
period of the planet was obtained. A few more observations were 
made in April and July 2014 to characterize the observed long¬ 
term trend. Simultaneous Th-Ar calibration was used to achieve 
high RV precision. The complete set of observations is shown 
in Table |T| When possible, we observed the fast-rotating, hot 
star 77 UMa (B3V, V=1.84, Vsin/' = 205 kms -1 ) immediately 
afterwards; this was our tool to remove telluric lines. 


Table 1 . HARPS-N RV observations of t Boo. This table is 
available in its entirety online at the CDS. 


BJD UTC -2450000. 

RV (ms ‘) 

RV err (m s 

*) Bis. span (ms *) 

6396.517465 

-16145.21 

0.83 

-104.71 

6396.518460 

-16136.13 

0.83 

-114.74 

6396.519444 

-16141.63 

0.84 

-109.08 

6396.520439 

-16141.12 

0.80 

-102.26 


To study the star-planet interaction, we needed to reach a 
high signal-to-noise ratio (S/N) near the Can H&K lines. To 
match this requirement with a single exposure would cause the 
saturation of the rest of the spectrum, given the spectral energy 
distribution of r Boo and the lower efficiency of HARPS-N in 
the blue region. The need to avoid saturation was combined with 
the requirements of an asteroseismic feasibility study of r Boo 
with HARPS-N, allowing us to make a synergy between differ¬ 
ent science themes of the GAPS campaign. Our observational 
strategy consisted of taking several one-minute exposures and 
then averaging them to obtain a single spectrum with a very high 
S/N. In this way, we obtained both the high-cadence spectra to 
monitor the solar-like oscillations and the mean high S/N spec¬ 
tra to study the Ca n H&K variations over the orbital period of 
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the planet. Moreover, we also obtained RV values to study the 
long-term variability of the system. 

We developed an interactive program written in Python 
dBorsa et al.ll2013l) to average the raw images pixel by pixel, ad¬ 
justing the header of the created files. The resulting mean FITS 
files are ready to be passed through the HARPS-N pipeline. In 
this way they are reduced exactly in the same way as all the FITS 
files acquired with HARPS-N (Table[2}. 

The Julian dates of the mean files calculated by the pipeline 
are not the correct dates because of the overhead time between 
the single exposures.This time is ~25 sec, therefore its influence 
on consecutive 1 min exposures is relevant. We corrected for this 
by taking an average of the Julian dates of the single images, 
weighted on their respective S/N: in this way, we introduced a 
sort of exposure-meter information in the mean file. The RVs 
were then corrected for the change in the barycentric Earth radial 
velocity (BERV) between the pipeline-estimated Julian date and 
the corrected date. For this purpose, we created a new tool and 
verified it to be co mpara ble with the HARPS-N pipeline at the 
level of ~ 3 cms -1 dBorsa et al.ll2013i) . 

Table 2. RVs for the r Boo nightly mean spectra. (f> refers 
to the orbital phase of the planetary companion, based on the 
ephemeris of Sect. 0 and considering the planetary inferior con¬ 
junction as 0 = 0 . 


Night 

BJD utc -2450000 

RV [ms" 1 ] 

RV err [ms 1 ] 

0 

1 

6396.531788 

-16135.08 

0.86 

0.73 

2 

6397.536059 

-16716.81 

0.77 

0.03 

3 

6398.499820 

-17016.77 

0.93 

0.33 

4 

6399.548081 

-16231.25 

1.39 

0.65 

5 

6401.498724 

-17085.13 

0.79 

0.23 

6 

6402.497328 

-16501.11 

1.31 

0.53 

7 

6406.653040 

-16169.98 

0.78 

0.79 

8 

6407.703825 

-16904.41 

0.89 

0.10 

9 

6408.686713 

-16878.51 

0.96 

0.40 

10 

6410.680969 

-16603.65 

1.34 

0.00 

11 

6421.517542 

-17054.54 

0.97 

0.27 


The data (324 single exposures and 11 mean exposures) 
were finally reduced using the HARPS-N data red uctio n soft- 
ware (DRS) pipeline on the Yabi platform. Yabi dHunter et al.l 
l2012l> is a Python web application installed at IAi0 in Trieste 
that allows authorized users to run the HARPS-N DRS pipeline 
on their own proprietary data with custom input parameters. 

The pipeline installed at the TNG estimates the radial ve¬ 
locities of the targ e ts by com puting a cross-correlation func¬ 
tion (CCF, iPepe et alTI 120021) using the best-suited line mask 
of the available masks (G2, K5, or M2 spectral type). Taking 
advantage of the Yabi platform, we were able to c reate and im¬ 
plement a new custom mask for r Boo (iRainerl [201 3t iGrattonl 
12013l l. Using the standard G2 mask as a starting point, we mea¬ 
sured the depths of several unblended lines in a well-exposed, 
high S/N r Boo spectrum. From these measurements we found 
the following empirical correlation between the line depths (LD) 
of the G2 mask and those of the r Boo spectrum: LD T boo = 
0.0409321 • e 2 - 806 54 t.OG 2 pj s i n g this equation, we corrected the 
values of the line depths in the standard G2 mask and created 
the new custom mask (3625 photospheric lines) that we used for 
the data reduction. In addition to this, we increased the width of 

3 http://ia2.oats.inaf.it/ 


the half-window of the weighted CCF from the default value of 
20 up to 30 km s (about twice the value of Vsin i of the star) to 
clearly cover the continuum around the wings of the mean line 
profile. These improvements reduced the RV errors by ~5%. 

3. Spectral analysis 

The analysis of the magnetic activity of r Boo was performed us¬ 
ing the nightly high S/N mean spectra. The CCF provided by the 
DRS pipeline was computed by cross-correlating the spectrum 
of t Boo with the custom mask (see Sect. 0 : it is thus essen¬ 
tially a high S/N mean line profile, and it is highly sensitive to 
the presence of active regions on the stellar photosphere. 

To study the variability of the photospheric activity, we mon¬ 
itored the line profile variations of the CCF, in particular its con¬ 
trast (CCFc) and full-width at half-maximum (FWHM), the bi- 
sector inverse slope (BIS, [Gravl 12008}) , and the v asy parameter 
dFigueira et al.l(2013h . We computed these indicators by means 
of custom-made routines, developed by our team for HARPS-N 
spectra. 

We also monitored the variability of the chromospheric ac¬ 
tivity. We focused on the Can H&K lines (3968.47 A and 
3933.66 A), the Na i D [2 doublet (5889.95 A and 5895.92 A), the 
He i D 3 triplet (blend at 5875.62 A), and the H„ line (6562.79 A) 
as chromospheric diagnostics. 

The chromospheric indicators were analyzed starting from 
the HARPS-N one-dimensional spectra reduced with the DRS 
pipeline, which computes the localization of spectral orders 
on the two-dimensional images, performs order extraction and 
stitching, corrects for flat-fielding, rejects cosmic-rays, and cal¬ 
ibrates the wavelength. Our analysis is based on the differen¬ 
tial comparison between the spectra in the series realigned in 
the wavelength space. This approach allows us to avoid the un¬ 
certainties related to the absolute flux calibration and contin¬ 
uum norm alization and is therefore m ore robust than the method 
adopted bv IScandariato et ahl (12013l) . which depends on the nor¬ 
malization to the continuum (a difficult task in the presence of 
broad lines such as the Ca n H&K and order stitching). 

Here we discuss how we reduced the data to perform the dif¬ 
ferential analysis, while in Sect. 16.11 we analyze the variability 
in the spectra of r Boo. In the following analysis we excluded 
the data point of night 4 because it has low S/N, and no spectro¬ 
scopic standard star was observed because of bad weather. 

The Na i D 12 doublet, the He i D ? triplet, and the H„ line are 
affected by telluric contamination, which is variable from night 
to night as a result of different airmasses and sky conditions. 
For each night of observation, telluric features were removed 
from the target spectrum by comparing it with the spectrum of 
the standard star 77 UMa. We corrected for telluric absorption 
using the task telluric in IRAF (Tody 1993). In contrast, the 
Ca 11 H&K spectral region of the standard star does not show any 
telluric feature above noise. Telluric correction was therefore not 
performed on the corresponding spectral interval of r Boo to 
preserve the original S/N. 

To carry out the spectral differential analysis, we needed to 
rescale all the spectra to the same flux scale. To this purpose, 
we selected the spectrum at mid-time (i.e., the spectrum of night 
6 ) with respect to which we performed the differential spectral 
analysis. This choice aims at minimizing any time-dependent in¬ 
strumental effect along the time series of spectra (see below). 

For each diagnostic we extracted the corresponding spectral 
region, and we compared each spectrum in the series with the 
reference spectrum on a pixel-by-pixel basis (Fig.|T). taking ad¬ 
vantage of the wavelength stability of HARPS-N. In this way. 
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Table 3. Spectral ranges used for the analyzed chromospheric 
diagnostics. 


Line 

Vacuum wavelength 
(A) 

Spectral width 
(A) 

Core width 
(A) 

CanK 

3933.66 

3.5 

0.5 

CanH 

3968.47 

3.5 

0.5 

Nai D 2 

5889.95 

3.5 

0.3 

Nai Dj 

5895.92 

3.5 

0.3 

He i D 3 

5875.72 

3.5 

0.6 

H„ 

6562.79 

3.5 

0.5 



Figure 1 . Example of the flux calibration of the mean spectrum of 
night 3 over the mean spectrum of night 6 in the Can K spectral range. 
The best fit (red line) is computed using the black dots, which represent 
instrumental fluxes out of the line core. The instrumental fluxes of the 
line core (in gray) are not included in the fit. The equation of the best fit 
is also shown. 


the reference spectrum defined the common flux scale, and all 
the other spectra were rescaled accordingly using a linear best fit 
of the flux vs flux relation. The intercept adjusts the background 
level, while the slope rescales intrinsic stellar fluxes. Line cores 
were excluded from the best fit to avoid any intrinsic nonlinear¬ 
ity that might have been introduced by chromospheric variabil¬ 
ity. The widths of the analyzed spectral ranges and the line cores 
are reported in Table [3] 

This procedure generally led to satisfactory results, that is, 
fluxes were well aligned along a straight line (Fig. Q}. Still, in 
a few cases especially in the Can K line, the residuals of the 
fit were higher by ~5% than the noise over small wavelength 
ranges, typically <1 A (Fig. [2}. These distortions, frequently ob¬ 
served during the first year of operation of HARPS-N, are not 
relevant in the measurement of RVs, but introduce spectral inho¬ 
mogeneities on a night-to-night basis that may hamper the SPMI 
analysis. To correct for them, we divided the spectra pixel by 
pixel by the reference spectrum, and we locally fitted the ratio 
with a low-order polynomial, excluding a narrow window cen¬ 
tered on the line cores where the SPMI signal might be mistaken 
as an instrumental effect. Finally, each spectrum was divided by 
the fitted ratio to remove the low-order flux variations. In Fig. [2] 
we show an example of the low-order correction in the Ca n K 
line. 

After removing the distortions, we averaged the time series 
of the corrected mean spectra for each spectral range. Then we 
computed the relative residuals for each mean spectrum with re- 



3930 3932 3934 3936 

UA) 


Figure 2. Example of the low-order correction of the Can K line. The 
black line is the ratio of one spectrum in the series to the reference 
spectrum, while the red line is the low-order fit (95% confidence band 
is shown with red dashes). The vertical dashed lines bracket the line 
core, in gray, which is excluded from the fit. 




Figure 3. Residuals relative to the average spectrum of the Can K (top 
panel) and H (bottom panel) lines. The color code (online version only) 
is reported in the legend, together with the S/N in the line for a 1 A 
spectral element (night 4 has been excluded, see text). The residuals 
have been smoothed to avoid cluttering; for this reason, they slightly 
exceed the vertical dashed lines. 


spect to the averaged spectrum and integrated them over a ~1A 
window around the line center to obtain the integrated relative 
deviation (IRD). Figure [3] shows the residuals, which have been 
smoothed to show the variability in the Ca n H&K line core more 
clearly. We computed the uncertainties on the IRDs by summing 


4 





















































F. Borsa et al.: Putting exoplanets in the stellar context: r Boo A 


in quadrature the photon noise and the uncertainties on flux cal¬ 
ibration and low-order correction. In this way, we take into ac¬ 
count the same sources of error more than once, thus obtaining 
an upper limit to the true uncertainty. 

We note that the relative residuals in the line cores are com¬ 
parable to or slightly smaller than our low-order corrections to 
the spectra. The reliability of our approach rests on the verifi¬ 
cation that the recovered Can H&K signals (IRDs) are reason¬ 
ably correlated, and the same holds for Ca n H+K vs H„ IRDs 
(Sect. ED. We have verified that our low-order correction im¬ 
proves these physically meaningful correlations. 

4. Stellar parameters 

We derived the physical parameters of r Boo by means of 
fittin g sy nth etic spectr a (ATLAS9 stellar atmosphere models, 
ICastelli & Kurucal2004lf t o a normalized HARPS-N mean file 
using the SME software bv IValenti & Piskunov I (119961) . The hy¬ 
drogen lines could not be used to determine the temperature be¬ 
cause they span several spectral orders, which introduces prob¬ 
lems in the normalization of the spectra. Therefore we used the 
region between 5160 and 5190 A, where the Mg triplet is found. 

We computed the Vsin i using the Fourier transform of the 
CCF of a mean spectrum (see Sect. 14.11) . This value was kept 
fixed in the fit while the other parameters were left free to 
vary. Table [4] shows our results, which agree well with the 
values found in the literature (e.g.. ISantos et al.ll2013l) . A por¬ 
tion of the observed and synthetic spectra is shown in Fig. [4] 
An independent check of the stellar parameters was made with 
the metho d based on the equivalent widths of s pectral absorp¬ 
tion lines (ISousa et al. li2007tlBiazzo et al. 1120121 and references 
therein). Even if this method is not preferable because of the 
high Vsin i value of t Boo, results were compatible with those 
estimated with the synthetic spectra method. 


Table 4. Stellar parameters derived for r Boo. 


Parameter 

Value 

Teff [K] 

6399 ± 45 

log CJ [cuts -1 ] 

4.27 ± 0.06 

[Fe/H] 

0.26 ± 0.03 

Vsin i [km s~*] 

14.27 ± 0.06 

Luminosity [L 0 ] 

3.06 + 0.16 

Mass [M e ] 

1.39 + 0.25 

Radius [R 0 ] 

1.42 + 0.08 

V breakup [km S ] 

352 ± 67 


We determined stellar radius and mass (Table [4} by com¬ 
bining the spe ct rosc opic results with the parallax (64.03 + 0.19 
mas, Ivan Leeuwen I 20071) and using a bolometric correction 
t-0.0078 7ibrres ll201(*lHowerlll996f) . 


4.1. Differential rotation 

We found evidence of solar-like differential rotation in r Boo 
from studying the Fourier transform of the mean line profiles of 
the HARPS-N mean spectra. We used two different mean line 
profiles: the CCF computed by the HARP S-N pip eline and the 
profile obtained using the FSD software ( Donati et al. I [ 19971 ) 
on the wavelength regions 4415-4805, 4915-5285, and 5365- 
6505 A. The CCF was computed using the Yabi interface, which 



Wavelength (angstroms) 


Figure 4. Fit of the observed mean spectrum (solid black line) with the 
synthetic spectrum (dashed red line) in the wavelength region consid¬ 
ered, using the SME software. 



Figure 5. Amplitude of the Fourier transform of a mean line profile. 
Red dashed line refers to the LSD software, black solid line to the 
HARPS-N CCF. The horizontal lines show the noise level. 


means that we were able to reduce the spectra with our custom 
mask, and we further increased the CCF half-window width up 
to 40 kms -1 , to ensure that we captured the whole mean line 
profile and the continuum. 

We used a Fourier transform on both the LSD and the CCF 
profiles (Fig. [5} and found the values of the first two zero po¬ 
sitions q\ and qi. The q\ position gives an estimate of the pro¬ 
jected rotational velocity Vsin/, while the ratio q?/qi is a n in¬ 
dicator of differential rotation (iReiners & Schmitt! [2002]) , ei¬ 
ther solar-like (q 2 /q\ < 1.72) or antisolar (q 2 /qi > 1.83). We 
found qilq\ — 1.53 + 0.08 for the LSD profile and q 2 /q\ — 
1.63 ± 0.04 for the CCF. The larger error in the LSD profile 
arises probably from the fact that we did not use the whole spec¬ 
tral range. The se values are compatible w ith the results from 
[Reiners Id® (q 2 /qi = 1.57 + 0.04) and lCatala et all d2007h 
(q 2 /qi = 1.60 + 0.02). 

The q 2 /q\ we found can indicate either solar-like differential 
rotation (equator rotating faster than the poles) or strong gravity- 
darkening in a rigidly rotatin g star. In the second case, the fol¬ 
lowing empirical equation ([Reiners 1120031) can be used to infer 
the rotational velocity of the star: 

q 2 /q\ - 1-75 + civ + bv 2 , (1) 
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where a and b are parameters that depend on the stellar spec¬ 
tral type. In particular, these parameters are a = 0.172 x 10 5 
and b = -0.993 x 10 ' 6 for F0 stars and a = 0.184 x 10 ' 3 and 
b = -0.116 x 10 ' 5 for GO stars, r Boo is classified as an F 6 
star, falling in the middle of these values, which means that to 
justify the values of q^/qx found by means of rigid rotation, the 
star’s rotational velocity should be either between 415 and 450 
km s ' 1 (in the CCF case) or between 520 and 550 km s ' 1 (in the 
LSD case). Considering the values of mass and radius of TableE 
we can compute the breakup velocity breakup = 352 km s' 1 . 

For the values qilq\ to be caused by gravity-darkening, the 
star would have to rotate faster than the breakup velocity. We can 
therefore reliably establish that the two qi/q\ values found show 
solar-like differential rotation in r Boo. 

Adopting an inclination angle for the star i = 40° 

(iDonati et al.120081 ) and using the qilq\ values, we computed the 
differential rotation parameter a - 7 ^,where f 2 () is the equato¬ 
rial angular velocity of the star and AO is the difference between 
the equa torial and polar angular v elocities. Using the empirical 
relation dReiners & Schmitt]l2003l) 


a 


2.74-5.16 



+ 4.321— 
<7i 



( 2 ) 


we found a = 0.24 ± 0.07 from qilqx — 1.53 (LSD profile) 
and a = 0.16 + 0.04 from < 72 / 1/1 = 1.63 (CCF) , both compatible 
with the value a = 0.18 found by ICatala et all (120071) . The LSD 
resul t is exactly the same as the differential rotation found by 
IDonati et al.l ( 2008 ) using spectropolarimetric measures. 

Our work shows that the LSD profile and the CCF provide 
similar values of qilqi- As such, the CCF computed internally 
by the HARPS-N pipeline can be used as an indicator of differ¬ 
ential rotation, at least for high S/N spectra. 


5. Orbital fit 

We used data from Table [T] to perform the orbital fit, but when 
multiple exposures were taken during one night, we used the RV 
value corresponding to the nightly mean spectra (e.g., Table[2} to 
be independent of stellar oscillations. The final HARPS-N set of 
measurements is composed of 20 RV points. In addition to our 
HARPS-N data, we used the recently released archival data of 
the Lick Observatory dFischer et al.ll2014» . These data were ob¬ 
tained with four different iodine-cell setups (identified as number 
2, 13, 6 , and 8 in Table [5} that were taken in consideration dur¬ 
ing the analysis adding a RV shift as a free parameter to each set. 
We excluded all the data with cr > 30 ms ' 1 (e.g., all the data of 
Lick setup number 2) and those that did not maintain the same 
instrumental setup continuously. We used 166 Lick RV measure¬ 
ments in the orbital fit. By combining Lick and HARPS-N data, 
we obtained a dataset composed of 186 measurements, spanning 
a time interval of about 20 years. 

We fitted the RVs w ith a Levenberg-Marquardt algorithm 
JWright & Howard] 120091) . using a two-planet model (Fig. Eb 
to take into account both the short-term variability due to the 
known planet and the long-term trend due to the stellar binary 
companion (Sect. 15.1b . To estim ate the error-bars, we used a 
bootstrapping code dWang et al.ll2012l) . The orbital parameters 
obtained are reported in Table 0 

We analyzed the RV residuals shown in the bottom panel 
of Fig. [ 6 ] in frequency, but did not find any trace of additional 
periodicity.We note in particular that RV residuals of HARPS- 
N show a clear correlation with the pipeline-estimated bisector 



- 0 2000 4000 6000 

: [r;:,lQ-?4500C0] 


Figure 6. Two-objects orbital fit (gray) for the Lick Observatory 
archival data (red upward triangles, pink squares, and green downward 
triangles corresponding to the different setup configurations 13, 6, and 8 
respectively) and HARPS-N data (blue circles). RV residuals are shown 
in the bottom plot. 


Table 5. Orbital parameters estimated for the r Boo system. Period and 
periastron of the binary companion r Boo B were fixed to those of the 
astrometric solution reported by Drummond (2014). 


Parameter 

Value 

r Boo b 

Period [days] 

3.3124568 ± 6.9 x 10' 6 

T periastron [BJDuTC-2450000] 

6400.94 ± 0.30 

K [ms' 1 ] 

471.73 ±2.97 

e 

0.011 ±0.006 

0 ) [deg] 

113.4 ±32.2 

y [ms' 1 ] 

0.0 

m sin i [Mj up ] 

4.32 ± 0.04 

semi-major axis [AU] 

0.049 ± 0.003 

r Boo B 

Period [years] 

964 (fixed) 

T periastron [BJD UTC -2450000] 

12670 (fixed) 

K [ms' 1 ] 

1217.06 ±222.36 

e 

0.71 ±0.22 

OJ [deg] 

94.1 ±64.0 

y [ms' 1 ] 

-1099.5 ±273.0 

m sin; [M G ] 

0.4 ±0.1 

semi-major axis [AU] 

109 ±7 

offset Lick,s e tupi3 [m s ] 

0.0 

offset L , c /^ erMp6 [m s ] 

-23.17 

offset L ; c ^/ w/ , 8 [m s ] 

19.03 

offset harps-n [ms' 1 ] 

-16294.87 

independent long term trends 

RV slope,_, rfvS(; ,„ p2 [ms' 1 y' 1 ] 

2.7 ±5.5 

RV slope,. Ict , M ., up i3 [ms' 1 y' 1 ] 

-13.8 ±5.6 

RV slope,irt.sempc, [ms' 1 y' 1 ] 

-8.8 ±3.4 

RV slope,[ms' 1 y' 1 ] 

-18.6 ±2.4 

RV slop e HA Rps-N [ms' 1 y' 1 ] 

-36.6 ± 4.0 


span (Fig. [7]), which supports their stellar activity nature (dis¬ 
cussed below in Sect. [ 6 }. 
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Figure 7. Residuals of the planetary fit for the HARPS-N measure¬ 
ments vs the pipeline-estimated bisector span. The green squares refer 
to exposures taken at JD 2456421, during which the target was observed 
between clouds. Red circles are the measurements for the nightly mean 
spectra. Error bars for the bisector span are taken as twice the RV error. 



5.1. Stellar companion 

Many astrometric measurements have been made to estimate the 
orbital parame ters of the binary stellar system, with th e most re¬ 
cent solutions dRoberts et al. l20lHiDrummondll2014ll claiming 
that the binary companion has a period < 1000 years and will 
approach periastron within the next two decades. 

To check whether the RV data match the astrometric best-fit 
orbit well, we treated each Lick instrumental setup and HARPS- 
N separately. We fitted a one-planet model, imposing that the 
orbital solution for r Boo b be fixed, and keeping a long-term 
linear trend free (lower part of Table 0. In this way, we were 
able to show the change of the slope in RVs caused by the binary 
companion r Boo B for the first time. The information of the 
astrometric_orbital parameters dDrummond 1 1 20 14 3 and parallax 
(Ivan Leeuwen 1 1200% indicates a mass sum of 1.8Mq and thus 
a value of 0.4Mq for r Boo B. We used this value to calculate 
the astrometric-based orbital RV slope caused in r Boo A by the 
stellar companion (Fig. [8]). 

Clearly, the data are insufficient for more constraints, but the 
best-fit astrometric solution agrees well with the present RVs. 
The fact that the slope calculated based on HARPS-N data is 
more than twice those calculated based on Lick data demon¬ 
strates that the star is rapidly accelerating as a result of the ap¬ 
proaching periastron. A monitoring of the RV trend during the 
next years will enable us to obtain more reliable orbital parame¬ 
ters for the stellar binary companion. 


6. Stellar activity 

In this section we aim to study the SPMI in the t Boo system. To 
this purpose, we extensively analyze and compare all the mag¬ 
netic activity indicators available, both computed directly from 
the spectra and extracted from the CCF (see Sect. [3]). 

6.1. Analysis of the IRD vs IRD relationships 

6.1.1. Can H&K and H tJ , 

Figure [3]shows that the cores of the Ca n H&K lines are affected 
by variability that is stronger than noise in the corresponding 
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Figure 8. Points represent the values for the long-term RV slopes fit¬ 
ted on Lick and HARPS-N data. The red square is the value for the first 
Lick setup with large error bars, which was not taken in consideration 
when we fit the planetary parameters. The black dotted-dashed line rep¬ 
resents the fitted orbit of the binary companion from the two-objects fit 
of Table[5] The blue dashed line is the long -term slop e for the best-fit or¬ 
bit for the stellar companion reported bv lDrummond I ( 120141) . assuming 
a mass value of 0.4 Mq. 



spectral regions. In Fig. [9] (left panel) we plot the correspond¬ 
ing IRDs of the two lines (IRDh and IRDk): we found that the 
two proxies are strongly correlated with each other with a confi¬ 
dence level of ~94%; this shows that the detected signal is not an 
artifact produced by our data reduction, but arises from the mag¬ 
netic activity of the star. The significance of the correlation was 
computed with 10000 random permutations of the data, which 
were also randomly perturbed by the corresponding uncertain¬ 
ties. Since the uncertainties on the variables are similar, the best- 
fit line was compu ted by me ans of a ranged major axis (RMA) 
regression dLegendre & Legendrefll983h . 

Since the two H and K proxies are strongly correlated, 
we averaged them to obtain the collective IRDhk- We found 
that the latter stron gly correlates with the canonical l o g /^'hk 
( N oves et al.l 1984 1 computed following iLovis et al.l (1201 ill 
(Fig .[9] middle panel). 

Moreover, as shown in Fig. [9] (right panel), IRDhk also 
strongl y co rrelates with IRDh ■ A similar result was obtained by 
iFares et akl d2009l) , who spectroscopically monitored the activity 
level of r Boo in 2008. This is consistent with different studies by 
other author s dMeunier & Del fossc 2009; Martmez^AmaizeTdJ 
I20TqL I20TT1 IStelzer et al.l 120131 Kiomes da Silva et al.l 12014 . 
who found strong pairwise correlations between the H&K and 
the H, r lines despite the different formation heights (lower and 
upper chromosphere) and adopting different approaches (e.g., 
single-epoch comparison of many stars, multi-epoch variability 
study of single objects, or emission line fluxes vs indices as rel¬ 
ative measures). 


6.1.2. Na 1 Di > 2 doublet and He 1 D 3 triplet 

We found that the IRDs of the Nai Di and Nai D 2 lines corre¬ 
late linearly, but with a confidence level not better than ~92% 
(Fig. [TO] left panel). We carefully checked the spectra and did 
not find any hint of a poor telluric correction or uncorrected ter¬ 
restrial emission lines above the noise level. As described above 
for the H&K lines, we averaged the two proxies to obtain the 
IRDdi 2 indicator. 
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Figure 9. Left panel: IRD H vs IRDk- Middle panel: IRDhk vs LoglRJ^). Right panel: IRDhk vs IRDh„ . In all panels, the color code (online 
version only) is the same as in Fig. [3] while the number at the bottom is the confidence returned by the correlation test. The numbers above the 
data points indicate the ordered sequence in the series of observations. Gray dashes show the RMA best-fit. 


We did not find any correlation from plotting IRDdi 2 against 
IRDhk (Fig . ITOl middle panel). This is consistent with the re¬ 
sults of iDfaz et al.l d2007l) . who only found a good correlation 
between the Na and Ca doublets for stars with Balmer lines in 
emission (not the case of r Boo). In any case, it is interesting to 
note that there is a bifurcation in the sample of measurements 
(Fig. [TOJ middle panel). 

In a similar way, the IRDdi2 vs IRDq3 plot seems to clus¬ 
ter into two subgroups (Fig. ITOl rightjpanel). The origin of the 
two trends is still unclear. iLandmanl (1198 lh . analyzing a few qui¬ 
escent prominences on the Sun, obtained a similar bifurcation 
by comparing the intensities of the Hei D 3 triplet with the Nai 
Di and D 2 lines. He conjectured, with no clear demonstration, 
that since the He triplet is back-heated by coronal UV radiation, 
then the two branches reflect the association of prominences 
with nearby coronal activity of different levels. The author re¬ 
marked that the origin of the two branches is related to the He 1 
D 3 triplet, as it shows the same bifurcation when compared to 
the intensity of the Ca 11 T8498 line, while the latter well corre¬ 
lates with the Na 1 D 12 lines. In contrast, we found the bifurcation 
when IRDdi 2 was compared with both IRDd 3 and IRDhk, while 
there is no bifurcation in the IRDd3 vs IRDhk plot. Thus our 
results indicate that the formation of the Nai doublet should be 
investigated to explain this twofold behavior. 

Another possibility is the presence of prominence-like struc¬ 
tures around the stars, formed with matter evaporating from the 
planet and supported by the magnetic field of the star ( Lanza! 
12014) . Since exoplanets are expected to be richer in metal than 
stars (e.g., iFortnev et al.l2006l) . this may explain why Na 1 shows 
the two branches instead of the He 1 D 3 triplet. 

Focusing on the He 1 D 3 triplet, we found that the scatter of 
the IRDd 3 measurements is slightly larger than the uncertain¬ 
ties, indicating that weak variability in this proxy has occurred, 
if any. Moreover, we did not find the correlation with the_Can 
H&K lines that has previously been found bv lGarcia-Lonez et alJ 
(119931) in their sample of F-type main-sequence stars. This is 
probably because their statistical sample spans a wide range of 
log R'hk, while r Boo’s chromospheric variability remained too 
limited during our observational campaign to detect such a cor¬ 
relation. 


6.1.3. Parameters of the CCF 

The very high S/N of the mean spectra allowed us to study the 
CCF variability with a high degree of accuracy. In particular, 
since the CCF is obtained from a large set of photospheric lines, 
the CCF variability is indicative of the source and level of pho¬ 
tospheric activity. It is thus interesting to study the CCF profile 
to understand the phenomenology of magnetism on the surface 
of r Boo and how it is related to chromospheric activity. 

As suggested by iNardetto et alJ ( 20061) . one way to analyze 
the CCF profile is to retrieve its width and contrast to the contin¬ 
uum by means of a Gaussian best-fit. They also introduced the 
possibility to fit asymmetric profiles by tweaking the Gaussian 
model. For r Boo, both methods led to poor fits of the observed 
CCF profiles, probably because the differential rotation com¬ 
bined with the rotation rate of the star broadens the line profile 
in a non-Gaussian way. 

We therefore approached the CCF profile analysis by averag¬ 
ing the CCFs returned by the DRS from the series of ten nightly 
averaged spectra, obtaining a master CCF profile with a S/N 
higher than each nightly averaged CCF’s S/N. Then, we fitted 
each nightly averaged CCF with the master CCF, allowing for 
changes in FWHM and contrast (CCFc). 

We als o comp uted the BIS and the v asy parameters defined 
bv lFigueira et al.l ( i2013l) . The former measures the velocity dif- 
ference betw een the midpoints at the top and bottom of the CCF 
( Gray 2008 1. while the latter gives an estimate of the asymme¬ 
try of the CCF profile. By definition, the computation of both 
parameters is model independent. 

Since lOueloz et al.l (1200 lh . variations in BIS are the 
paradigm for the signature of activity-induced RV variations. 
These authors found a strong anticorrelation between the mea¬ 
sured RVs and the corresponding BIS of the CCF. Since then, 
the anticorrelation between BIS and residual radial velocities 
(RV„, v , the RVs minus the planetary best-fit orbit) is considered 
a clear signature of short-term (i.e., rotationally induced) photo- 
spheric variability. In our case, the BIS shows a high degree of 
anticorrelation with the residuals RV, ei (Fig. [7|i. This suggests 
that RV,,,< are due to stellar activity. 

Some other parameters show a strong anticorrelation: CCFc 
and FWHM (Fig. |TT] left panel). This is consistent with the fact 
that while the FWHM measures the width of the CCF at half 
maximum, CCFc counteracts this to preserve the area of the 
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Figure 10. Left panel: IRD D| vs IRD d , . Middle panel: IRDdi 2 vs IRDhk- Right panel: IRDdi 2 vs IRDdj. In all panels, the color code (online 
version only) is the same as in Fig.[3]and represents different observation dates. The numbers above the data points indicate the ordered sequence 
in the series of observations. In the left panel, the number at the bottom is the confidence returned by the correlation test. In the middle and right 
panels, the two dashed gray lines are drawn by hand to highlight the bifurcation discussed in the text. 
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Figure 11. Left panel: CCFc vs FWHM. Middle panel: CCFc vs BIS. Right panel: IRD d[ 2 vs CCFc. Symbols and colors (online version only) 
are the same as in Fig. [3] 


CCF. The two indices are thus sensitive to the broadening of the 
wings of the CCF. We infer, accordingly, that correlations be¬ 
tween these indices suggest a broadening variability of the CCF. 

Our results therefore indicate that the BIS follows the defor¬ 
mations of the CCF better than CCFc and FWHM which, in turn, 
are sensitive to the overall width of the CCF. The BIS and CCFc 
can then be considered the best representatives of two different 
families of indicators. Given the operational definitions of the 
BIS and CCFc, we can regard the former as an indicator of the 
mean distortion of the line profile and the latter as an indicator 
of the mean strength of the photospheric lines. 

The weak variability of BIS with respect to measurement un¬ 
certainties compared to CCFc (Fig. [TT] middle panel) may give 
some clues on the geometry of active regions on the stellar sur¬ 
face. The stability of the BIS and the variability of CCFc may in¬ 
dicate that the CCF changes its contrast with respect to the con¬ 
tinuum while preserving its overall asymmetry. This is consis¬ 
tent with a scenario with an active region around the pole of the 
star, which is not Doppler-shifted by stellar rotation and rapidly 
evolves in terms of contrast to the quiet photosphere (either in 
temperature difference and/or coverage factor). 

We did not find any significant correlation between v asy and 
the other parameters of the CCF. 


When we compared the CCF parameters with the chro¬ 
mospheric indicators discussed in this section, we found the 
strongest correlation between IRDqi 2 and CCFc (p=92.5%, right 
panel in Fig. HB. The same applies for IRDhk and IRDh„ with 
a lower confidence level (p<90%). We infer that the variability 
of the Nal D 12 doublet arises from changes in the line formation 
physics at the photospheric level, rather than from genuine chro¬ 
mospheric emissi on. The doublet is formed in the lower chromo¬ 
sphere dTripicchio et al.ll 199 71 and is thus more closely related 
to the magnetic activity of the lowest chromosphere and photo¬ 
sphere. Moreover, as stated above, the Na I Dn doublet is a good 
proxy for chromospheric activity only for mid-to-late type stars 
with emission Balmer lines, which is not the case of r Boo. 

6.2. Time series analysis 

Using the ephemeris of r Boo b reported in Table [3] we searched 
for the eventual phasing of each of the indicators discussed so far 
with the orbital motion of the planet. The most convincing cases 
are the phase-folding of IRDhk and IRDn a shown in Fig.[T2l(left 
and middle panel), where we set the phase cf> — 0 at the planetary 
inferior conjunction. In these diagrams we overplot the weighted 
least-squares best fit of the form IRD=A- sin 2n(f)+B- cos Itttft+C, 
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Figure 12. Left panel: Phase-folded diagram of IRDhk- Middle panel: IRDh„. Right panel: Time series of the CCFc measurements. The planetary 
inferior conjunction is at 0=0. Colors (online version only) and numbers are the same as in the previous figures. In the left and middle panel, the 
gray dashed line is the least-squares best fit discussed in the text. We also report in the text the relative likelihood of the sinusoidal fit with respect 
to the flat model returned by the AICc. 


where <p is the orbital phase and A, B, C are the parameters to 
be fitted. We excluded JD 2456396 (i.e., night 1) from the fit be¬ 
cause it clearly is an outlier. The coefficient of determination R 2 
of both fits is ~50%, which means that the model explains ~50% 
of the total variance of the sample. Moreover, the likelihood ratio 
test suggests that the sinusoidal fit is to be preferred with respect 
to the flat model with a relative likelihood close to 100%. 

For all the other indicators, both chromospheric and CCF- 
related, we obtained a lower R 2 coefficient, that is, we found 
weaker or no evidence of such a phasing. Conversely, we found 
that CCFc, and consequently IRDdi 2 , show a clear trend with 
time (right panel in Fig. IT2l >. If we interpret that CCFc anticor- 
rel ates w ith mag n etic a ctivity (either spot-like or plage-like, see 
iDumusaue et al.ll2014l) . then the decrease of CCFc between JD 
2456402 and 2456407 suggests that r Boo increased its activity 
level, and conversely it turned over toward a quieter configura¬ 
tion between JD 2456411 and 2456422. 


6.3. Discussion 

In Sect. 16.11 we argued that the CCFc shows variations on 
time scales longer than the stellar rotation. Polar magnetic re- 
gions on t Boo have already been detected in several epochs 
by iFares et all d2013l) . who remarked that the polarity of the 
magnetic field switches about every two years. Moreover, inside 
each activity cycle they found signatures of a rapid evolution of r 
Boo’s magnetic field. This supports the hypothesis that the vari¬ 
ations of CCFc are driven by a large-scale evolution of magnetic 
activity, even along time spans as short as ~20 days (Fig. El 
right panel). 

Previous photometric monitorings of the star have never 
detected photometric variability abov e the mmag level 
dBaliunas et al.l 1 1997 :: Walker et al.112008 ). Accordingly, we as¬ 
sume that the typical peak-to-peak variability is on the order of 
1 mmag (upper limit). With this assumption, and assuming that 
the active region is spot-dominated, we obtain that the active re¬ 
gion covers -0 ,1 % of th e vis ible hemisphere using the SOAP2.0 
code dPumusaue et alil2014i) . The corresponding peak-to-valley 
variations of RV, BIS, and FWHM returned by the simulations 
are <25 ms -1 , while those of our measurements are actually 
larger (~ 50 m s -1 , ~70 ms -1 , and ~150 m s _I respectively). 


Conversely, if we assume that the active region consists of 
a bright plage, then the coverage factor is ^2.5%. With this 
coverage, the simulated activity-induced signal in RV, BIS, and 
FWHM is on the order of 100 ms -1 , comparable with those of 
our measurements. We thus conclude that the high-latitude ac¬ 
tive region probably is plage-dominated. 

In Sect. 16.21 we briefly discussed the time variability of the 
most representative indicators, finding that the genuine chromo¬ 
spheric indicators, such as IRDhk and IRDh 0 , seem to be phased 
with the orbit al moti on of the planet. This has previously been 
found bv lWalker et al.l ( 2008 ). who detected a plage on r Boo at 
planetary phase 0 = 0.8 from Can H&K spectra in 2001 to 2003 
and claimed a photospheric spot using MOST photometric data 
taken in 2004 and 2005. The authors stated that the persistence 
of the active region at the same longitude is a strong indication 
that it is caused by a magnetic link between the planet and the 
star. On the other hand, Mathur et all (20141) have recently re¬ 
ported a few cases of bona fide single F-type stars observed with 
Kepler with active longitudes persisting on the stellar surface for 
many stellar rotations. 

Since r Boo’s rotational period equals the planetary period 
(with some degree of differential rotation, IFares et al.l2009t) . it is 
still unclear whether the phasing of chromospheric activity with 
the planet is due to SPMI or to an active region that simply coro¬ 
tates. In any case, if the SPMI scenario is confirmed, we remark 
that in our observing season the plage was located at <p ~0.1- 
0.2 (Fig. EK indicating that the magnetic connection between 
planet and star, if present, has moved between years 2005 and 
2013. This may be due to a poloidal reversal of r Boo’s mag¬ 
netic fields between these epochs. However, we remark that at 
phase 0.8 we found the largest scatter in the activity diagnostics: 
we have only two observations, and no statistical significance of 
this hypothesis can be assessed. 


7. Asteroseismology 

We have shown in the previous sections how the mean spectra 
can be used to derive very precise radial velocity measurements 
to study the planetary and binary orbits and very accurate indi¬ 
cators to monitor the stellar activity and star-planet interaction. 
We now describe the results obtained from the analysis of the 
high-cadence, short-exposure HARPS-N spectra. 
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Figure 13. Asteroseismology of r Boo: a linear fit was subtracted from the original HARPS-N data after the detrending from the planetary and 
binary orbits. 


We started our asteroseismic analysis with correcting the 
Doppler shift that is due to the orbital motion of the planet 
(Sect. [ 5 }. Then, to remove the residual low-frequency term left in 
the data, several approaches were tried, and we report the most 
powerful approach here: the subtraction of a linear trend from 
the data of each night. This procedure ensured that no spuri¬ 
ous low-frequency term (due to instrumental problems, plane¬ 
tary and binary orbits, or RV variations induced by activity phe¬ 
nomena with timescales longer than the duration of one observ¬ 
ing sequence) affects the radial velocity measurements that were 
used for the asteroseismic analysis. 

This procedure returned time series centered on 
RV=0.0 ms -1 on each night (Fig. IT3l >. The resulting RV 
curves show sporadic peak-to-peak variations of up to 15 m s _1 
(JD 2456399 and 2456410), but r Boo typically is much quieter, 
at lower than 10 m s -1 . This amplitude is larger than that of the 
solar twin 18 Sco dBazot et alj 2012 ), but comparable with that 
of the subgiant f3 Hvi ([Bedding et ail2007t) . 


The time series were then analyzed in frequency to detect 
short-scale periodicities. We used the iterative sine-wave least- 
squares method dVanicekll97lh and checked the results_with the 
generalized Lomb-Scargle periodogram (iZechmeister & Kursterl 
12009 1. The power spectra are very similar and strongly affected 
by the sampling of the observations. The spectral window shows 
strong aliases at multiples of 1 d _1 , i.e., 11.57 /rHz (Fig.[14j inset 
in the top panel). 


As expected because of the limited time coverage, the fre¬ 
quency analysis of the RV values cannot supply a unique de¬ 
termination of the asteroseismic content. The power spectrum 
is very noisy: a first relevant pattern occurs between 1.5 and 
1.8 mHz, followed by another between 2.4 and 2.9 mHz (middle 
panel). Both are affected by the spectral window effects: each 
pulsational mode of r Boo originates a structure similar to the 
spectral window, destroying the expected comb structure of the 
excited modes. The former pattern stands out a little more clearly 
over the noise, and it shows the highest peaks around 1.68 mHz 
(bottom panel). The corresponding period of 9.9 min can be 
glimpsed in Fig.|T3] it is traced by six consecutive measurements 
(60 sec exposure time, 25-35 sec overhead), but of course often 
modified by the interference with the other modes. The ampli¬ 
tudes of the peaks are very small, around 1.1 ms" 1 . This value 
takes into account the intrinsic incoherence (amplitude damp¬ 
ing, mode lifetimes) of the solar-like oscillations. We also cal¬ 
culated the level of the noise in the RV time series and obtained 
0.35 ms -1 . The corresponding S/N=3.1 leaves some uncertain¬ 
ties on the significance of the 1.5-1.8 mHz pattern since the 
threshold S/N=4.0 (e.g., [Bedding et al.ll2007l) was not reached. 
All these observational uncertainties are due to the limited time 
coverage: twice as many measurements would have allowed us 
to reach a threshold of S/N=4.0. 
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Frequency [mHz] 


Figure 14. Asteroseismology of r Boo. Top panel, spectral window 
of the HARPS-N data in the 0.0-4.5 mHz range and zoom in the 
0.0-0.5 mHz range (inset). Middle panel, power spectrum in the 1.0- 
4.0 mHz range. Bottom panel, power spectrum in 1.3-2.0 mHz range 
and power spectrum of this spectrum (inset). 


7.1. Asteroseismic results: observation vs theory 

The availability of a small set of high-precision radial velocity 
measurements allowed us to obtain an estimate of the asteroseis¬ 
mic parameters of r Boo making a very limited observational 
investment. We verify these results in a theoretical context. 

We computed the expected place of an excess of power in a 
star like r Boo. To do this, we used the scaling relations between 
stellar parameters (M, T e ff, and L\ Table 1) and the frequency 
of maximum power of the oscillations v max (e.g., IStello et all 
120091) . We obtained 1.98+0.46 mHz, compatible at +1 cr level 
with the observed value. We performed another check on the de¬ 
tection of v max by computing the power spectrum of the power 
spectrum to identify regularities in the detected frequencies. The 
most relevant feature was the 11.57 /(Hz spacing that is caused 
by the aliasing effect. After this, a peak at 94.4 /.(Hz has ap¬ 
peared (inset in the bottom panel of Fig. ITTb . If the structure 
centered at 1.68 mHz is due to solar-like oscillations, then the 
94.4 /(Hz spacing should be the large separation. The pair (v max , 
Av )-(168 0, 94) /( Hz matches the observed relation (see Fig. 2 in 
IStello et akll2009l) . Although not yet decisive, these results give 
us more confidence in the asteroseismic approach to study the 
stars hosting exoplanets by means of radial 

velocity measurements performed mainly for other purposes. 


8. Evolutionary stage 

It has been well demonstrated (e.g.. ICunha et al.ll2007 l) that an 
accurate determination of mass and radius combined with the 
spectroscopic measurements of atmospheric parameters allows 
us to clearly constrain the stellar structure of observed stars. 

To assess the information available from HARPS-N obser¬ 
vations, we therefore calculated a grid of theoretical structure 


models of r Boo evolved from a chemically uniform model on 
the zero-age main seq uen ce, u sing the ASTEC evolution code 
(iChristensen-Dalsgaardl[20b8al) and varying the mass and the 
composition to match the available atmospheric parameters (see 
Table H. 

The evolutionary models have been produced by employ¬ 
ing currenHphysical info rmatio n following the procedure de¬ 
scribed in lDi Mauro et al.l (1201 ll) . The input physics for the evo- 
lution calculations include d the OPAL 2005 equation of state 
('Rogers & Navvonovl2002l) . OPAL opacities (tlglesias & Rogers 
1996), and the NACRE nuclear reaction rates (lAngulo et ah 
19991) . Convectio n was treated accord ing to the mixing-length 
formalism (MLT; lBohm-Vitensel[T9 58) and defined through the 
parameter a - t/H p , where H p is the pressure scale height. 
The initial heavy-element mass fraction Z, was calculated from 
the iron abundance given in Table [4] using the relation [Fe/H] = 
log(Z/A) - log(Z/A) 0 , where (Z/A) is the value at the stellar sur¬ 
face. We obtained Z/A — 0,045 + 0.003 ass uming the solar value 
(Z/A) 0 = 0.0245 (iGrevesse & Noelsll 19931) . 

The resulting evolutionary tracks are characterized by the 
input stellar mass M, the initial chemical composition, and a 
mixing-length parameter. 

The location of the star in the H-R diagram identifies r Boo 
as being at the beginning of the main-sequence phase of core¬ 
hydrogen burning, with an internal content_of hydrogen in the 
core of about X c = 0.5. As predicted bv lLares et al.l d2009l) . it has 
a shallow convective region, with a depth of about D cz - 0.1/?, 
typical of stars in this phase. According to the stellar evolu¬ 
tion constraints, given the match with the observed atmospheric 
properties, and with the use of all the possible values of mass 
and metallicity, our computations show that the age of r Boo is 
0.9 ± 0.5 Gyr and the mass is M - 1.38 + 0.05M e . It is clear 
that the accuracy on the values of age and mass here inferred are 
dependent on the stellar model calculation procedure, and in par¬ 
ticular on the physical and chemical inputs (e.g.. lLebreton et ak! 
l2014l l. Here we simply note that lower values of the solar heavy- 
element abundance (Z/A) 0 led to a mass lower by 5% and an age 
higher by 30%. These values agree with those deduced above 
within 2-cr uncertainty limits. The lithium line at 6708A could 
not be used to confirm the young age, since T ? jj = 6399 K 
(Table[4} puts r Boo in the lithium dip (e.g jBalachandranl 19951) . 

The value of age obtained by direct modeling is better con¬ 
strained than the wide range of previous values obtained with 



To reproduce the oscillations obtained in Sect. [7] we se¬ 
lected the models that best fit the observations among all 
the computed models. Lor the selected models we calculated 
the a diabatic osc illatio n frequ encies using the ADIPLS code 
(Christensen-Dalsgaardi'2008h ). According to our calculations, 
this star should show solar-like pulsations with a spectrum in the 
interval 700 - 2600 /(Hz and frequencies equally spaced by a 
large separation of about Av = 95 ± 1 /(Hz. The observed large 
separation appears to agree excellently well with the theoretical 
separation. The range of frequencies is very similar, although the 
spectral window effects make a straight comparison very diffi¬ 
cult. We conclude that there is no discrepancy between the val¬ 
ues we inferred from the limited HARPS-N observations and the 
current modeling of solar-like oscillations in r Boo. 
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8.1. Tidal evolution 

lLanzal (1201 Ob and iDamiani & Lanza! (1201 5l) found that the ro¬ 
tational period P lot in a sample of main-sequence stars with 
7i.ii <; 6300 K accompanied by hot Jupiters with orbital period 
PQYb venfies the relationship 1 ^ Pvot/Porb > 2 (cf. Fig. 10 in 
iDamiani & Lanzal2015l) . Given its almost synchronous rotation, 
r Boo satisfies this relationship as well. The timescale of tidal 
synchronization of the stellar rotation is at le ast 1 - 2_orders_of 
magnitudes longer than the stellar age (e.g., iDonati et al]l2008t 
[Barker & Qgilviell2009 1t. suggesting that the system reached the 
ZAMS close to its current synchronous state. The subsequent 
evolution was characterized by a weak tidal interaction between 
the planet and the star with a very slow braking of the stellar ro¬ 
tation, owing to the lo w braking efficiency of th e stellar winds 
of mid-F type stars (cf. IDamiani & Lanzall2015l Sect. 6.1.2 for 
details). In other words, the present orbital and rotational angu¬ 
lar momenta in the r Boo system are probably close to the initial 
momenta with very little evolution since the system arrived on 
the ZAMS. 


9. Conclusions 

We studied the r Boo system by means of a new observational 
strategy applied to HARPS-N spectra that were collected in the 
framework of the GAPS project. This strategy allowed us to 
obtain both high-cadence one-minute exposures to monitor the 
solar-like oscillations, and high S/N unsaturated spectra to study 
stellar activity and star-planet interaction. 

We developed new computational tools to perform the co¬ 
addition of consecutive spectra and the subsequent average of 
them, preserving the correct computation of the BERV (and thus 
RV) values that have to be referred to the flux-balanced refer¬ 
ence time. Implementing a custom mask for r Boo allowed us to 
extract more accurate RV values; the adaptation of the Yabi plat¬ 
form gave us the possibility to test and verify these new tools in 
a user-friendly system. 

Our results are summarized in the following points. 

- We updated the ephemeris for the planet r Boo b and showed 
by means of our RVs that the binary companion r Boo B is 
rapidly accelerating, confirming the astrometric predictions 
that it is approaching the periastron on a highly eccentric or¬ 
bit. We plan to take additional observations in the next years 
to have a more reliable value of its eccentricity. 

- The SME software was applied to determine new stellar pa¬ 
rameters, giving the same results as the equivalent width 
method. In particular, we refined the values of r Boo’s mass 
and radius to 1.39 + 0.25 M 0 and 1.42 + 0.08 R 0 . We estab¬ 
lished that the star shows strong differential rotation, which 
we determined by means of both the LSD mean line pro¬ 
file and the CCF. We stress that the CCF computed by the 
HARPS-N pipeline can therefore be reliably used as an indi¬ 
cator of differential rotation. 

- The analysis of the correlations between several indices 
pointed out evident chromospheric activity. In particular, the 
activity indicators extracted from HARPS-N spectra suggest 
the presence of a plage around one of the poles of the star. 
The nature of the chromospheric activity remains uncertain. 
It is unclear if it is due to SPMI or to a corotating active 
region, or both. 

- Solar-like oscillations are detected in the RV time series. 
Although very limited by the spectral window, we inferred 
observational values of v max and Av. These values agree well 


both with the scaling relations and the asteroseismic model 
computed from our stellar parameters. This result supports 
our confidence in the application of the asteroseismic ap¬ 
proach to other bright stars hosting exoplanets to constrain 
their ages and masses. 

- From the evolutionary point of view, r Boo is at the begin¬ 
ning of the main-sequence phase of core-hydrogen burning, 
with an age of 0.9 + 0.5 Gyr. The model we built allowed us 
to further constrain the value of the stell ar mass to 1.38±0.05 
M 0 and thus, using i - 44.5 ± 1.5° dBrogi et al.ll2012l) . the 
mass of the planet to 6.13 + 0.17 Mj up . 
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